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}1^ ABSTRACT 

bJO 



We investigate the thermodynamical and causal consistency of cosmological models of the 
cubic Quasi- Topological Gravity (QTG) in four dimensions, as well as their phenomenological 

OO consequences. Specific restrictions on the maximal values of the matter densities are derived by 

requiring the apparent horizon's entropy to be a non-negative, non-decreasing function of time. 
The QTG counterpart of the Einstein-Hilbert (EH) gravity model of linear equation of state 
is studied in detail. An important feature of this particular QTG cosmological model is the 

fT) new early-time acceleration period of the evolution of the Universe, together with the standard 

late-time acceleration present in the original EH model. The QTG correction to the causal 
diamond's volume is also calculated. 
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1 Introduction 

Regions containing extremely dense matter are a common feature of all the big-bang inflationary 
cosmological space-times. The consistent description of such high energy states in the evolution 
of the Universe requires certain "higher curvature" extensions of Einstein-Hilbert (EH) gravity, 
involving powers (of the traces) of the Riemann tensor, believed to take into account the short 
distance quantum effects. Such corrections to the EH action are known to arise as counter- 
terms in the perturbative quantization both of matter in curved spaces [1] and of pure EH 
gravity as well [2,3]. The main problem with "higher curvature" gravity theories concerns the 
presence of higher derivatives of the metric in the equations of motion, which in general lead to 
causal and unitarity inconsistencies [2,3] when considered out of the framework of superstring 
theory. Nevertheless, the extensive studies of such models, in particular the so called modified 
/(i?)-gravities, have found interesting applications in different areas of modern cosmology (see, 
e.g., [4-6] and references therein). 

The present paper is devoted to the investigation of the effects caused by the higher curva- 
ture terms in a particular cosmological model in four dimensions, based on the simplest "most 
physical" extended gravity, given by the following cubic action for Quasi- Topological Gravity [7] : 

^GBL = /^^^^{i? - XL^ [2R^ - m^pR^P + 2R^p^,R'-P^^-] 
^3 + im^MSR^^^^Ri^^^^ _ mR^P^^R^^p^R^-^^ 

36i?"^'^ Rd'jRpS + ^R Rp^R'^ a\ +-Cmattcr}, (1) 
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where x^ = IGttG = 4/pj; l-p\ is the Planck length; A, ^ are dimensionless "gravitational coupling" 
constants, while L is a new length scale, which can be chosen as L = l-p\ by an appropriate 
redefinition of A and /i. We include the cosmological constant Aq > implicitly in the matter 
Lagrangian ^matter- 

The remarkable feature of this cubic extension of the EH action, is that the corresponding 
equations of motion for all conformally flat metrics, as for example those of domain walls and 
of flat Priedmann-Robertson- Walker (FRW) space-times: 

ds2 = -dt2 + a2(t)dxMxi , (2) 

are of second order [8] , while arbitrary metrics yield, in general, equations of motion of fourth 
order^. This fact, together with the introduction of an appropriate superpotential and the 



''The particular combination of quadratic terms represents the d — 4 Gauss-Bonnet topological invariant, and 
does not contribute to the dynamics. 



related BPS-like first order system of equations for the cubic QTG-matter model (1), derived in 
Ref. [8] , provide an efficient method for the analytic construction of a large family of exact flat 
FRW's solutions, representing asymptotically dS4 space-times. 

The most important and universal new property of these higher curvature cosmological 
models is that, differently from the EH case, the entropy prescribed to the apparent horizons [9] 

s{t) = ^^{l-2XL^H^ + 3fiL*H^) (3) 

with H = aja denoting the Hubble factor, is not automatically positive definite and increasing. 
As a consequence, the requirement of the thermodynamical consistency of these models: s{t) > 
and ds/dt > 0, introduces certain restrictions on the available maximal values ^max of the matter 
densities g = \H'^ < Qmax and certain minimal scales Lmin (related to ^max) up to which we 
can have a physically meaningful description of the Universe evolution within the framework of 
the cubic QTG cosmologies. 

In order to exemplify the effects caused by the higher curvature terms, we choose a par- 
ticularly simple and yet quite rich cosmological model, representing a QTG extension of the 
following EH model of a linear equation of state: 

Po/qo = wi- W2/go ■ (4) 

The matter content is that of a barotropic fluid with constant equation of state parameter wi, 
together with a dark energy "fluid" representing the cosmological constant Aq = x^W2/2(l+wi). 
This EH cosmological model has been widely studied [10-13], with a particularly remarkable 
result by Bousso et al. [10], who deduce the value of the cosmological constant for a universe 
dominated by dust through most of its history (wi = 0). The new QTG features established in 
the present paper are: The presence of a new (early-time) acceleration period; changes in the 
duration of the acceleration and deceleration periods, as well as of the future and past event 
horizon radii; and finally certain very small corrections to the volume of the causal diamond (to 
be compared with the EH one [10]). 

2 Modified FRW Cosmology 

Consider an universe filled with a barotropic fluid with energy density ^o and pressure po, 
components of a 'bare' energy-momentum matter-tensor TJiJ . For the ansatz (2), the QTG 
equations of motion derived from (1) are the modified Friedmann equations: 

x^^o = QH"^ (1 - f^L'^ H^) , (5a) 

x\qo + Po) = -4i? (1 - 3^L4 H^) ; (5b) 

^0 + 3H{po + qq) = . (5c) 

They reduce to the usual EH- Friedmann equations when /x = 0; otherwise the contributions from 
the QTG terms may be regarded as separating a "gravitational energy momentum tensor" T^ 
at the right-hand side of the Einstein equations, thus composing an effective energy momentum 
tensor Tf^, viz. 

r = T'^^ • T''^ = t(0) -I- t'^t^'^ 



where G^^ is the Einstein tensor. Notice that the effective energy-momentum tensor has all 

peff 



the properties of a perfect fluid tensor, i.e. its components, Tf^ = Diag {—Q,p,p,p) satisfy the 



usual, EH Priedmann equations: 

Q=^H\ g + p=-^H, (6) 

as well as the continuity equation, 

g + 3H{p + g) = 0, (7) 

which is a simple consequence of the Bianchi identities. Because of its direct connection to the 
Hubble function, the cosmological observations (of distances and red shifts) should perceive not 
the bare energy density qq, but rather the effective one, g, related to the former via Eq.(5a): 

go = e{i-fj.L^ig^) . (8) 

Although such a hydrodynamical interpretation of T^J; is rather formal, we further assume that 
this "effective fluid" obeys the weak energy condition, i.e. p + g > 0, what assures that H < 0. 
Then if the bare fluid also satisfies the weak energy condition, the function 

Co = 1 - 3^L4 H'^ = l- ^fiL^x'^g'^ (9) 

must be positive (or vanishing), as can be seen from Eq.(5b). While this condition holds auto- 
matically for ^ < 0, for positive values of fi the sign of Cq will depend on the value of the energy 
density; it vanishes for g = ■jj^\/i[\x and it is indeed negative for greater values of g. Since 
near a singularity the value of g grows without bounds, for ^u > it is only possible that both 
the bare and the effective fluids satisfy the weak energy condition in a nonsingular universe. 
This would be the case, for example, in bounce-like models beginning and ending at de Sitter 
spaces with (asymptotic) densities ^ds < it;^ vW/^- We leave the study of such spaces for a 
more thorough discussion in [14], and focus for the remainder of this letter in the singular cases, 
thus considering only /x < 0. 

3 Horizon entropy 

Killing event horizons are known to posses thermodynamical properties: a temperature T related 
to the surface gravity, and an entropy s which in EH gravity is given by the Bekenstein-Hawking 
formula s = A/ 4 (in geometrized units), A being the area of the horizon [15-17]. Then the 
Einstein equations can be rewritten as a Clausius relation [18], dE = T ds, for the flux of energy 
dE across the horizon^. 

The lack of time-like Killing vectors in non-stationary space-times - as for example typical 
FRW spaces - represents an obstacle in the definition of a surface gravity for the correspond- 
ing dynamical apparent horizons. A possible consistent generalization of the "Thermodynam- 
ics/Gravity" correspondence can nevertheless be achieved by using the Kodama vector to define 
the horizon's Kodama-Hayward temperature Tkh [25]. As it was recently shown by Cai et 
al. [21,22], the corresponding Friedmann equations, for EH gravity and for certain modified the- 
ories as well, turn out to be again equivalent to the Clausius relation dii^MS = ^KndsKW) with 
£^MS being the Misner-Sharp energy and skw an appropriately defined Kodama- Wald entropy, 
which for the cubic QTG cosmologies is given by Eq.(3).^ Notice that for de Sitter space-times. 



^Remarkably, this equivalence remains valid in a variety of "higher curvature" gravitational theories [19-22] 
with the horizon's entropy given then by the Wald formula [23, 24] . 

®The proof of the equivalence between the QTG modified Friedmann equations (5) and the above generalization 
of the Clausius relation is given in our forthcoming paper [26] . 



when H is constant and the apparent horizon coincides with the Kilhng event horizon, our 
formula (3) reduces to the known static QTG Wald entropy [9]. 

The consistent interpretation of sqtg(*)) given by Eq.(3), as an entropy for the apparent 
horizon in the considered QTG cosmological models requires that it must be a non-negative and 
non-decreasing function of time. The later is always true if both effective and bare fluids do 
satisfy the weak energy condition. Then as a consequence we have that s = — on i "a £• > 0. 
The restrictions imposed by the positivity condition sqtgC^) — ^^^ slightly more involved: 
depending on the signs and values of A and /i, they turn out to introduce certain upper bounds 
on the energy density g. 

Gauss-Bonnet Gravity: The topological nature of the GB term in the action renders the 
dynamics (i.e. the equations of motion) of GB gravity, for which fj, = and A 7^ 0, identical to 
that of EH gravity. But the horizon entropy is not simply proportional to H^^, in fact 

sit) = ^ {^ - 2XL^) . 

Thus if A < the entropy density is always positive, but for A > there is a value of H{t) - 
or, equivalently, of ^ - past which the entropy becomes negative. Therefore the assumption of 
positivity of entropy places as an upper boundary, qgb, on the possible values of the energy 
density: 

3 , , 

0< Q< QGB; gGB= ^212X ^'^^^ 

Quasi- Topological Gravity: Similar considerations applied to the QTG model, for the nega- 
tive values of the coupling ^ < we are interested in, lead us to the conclusion that the apparent 
horizon entropy is only positive for densities within the finite interval: 

< e < ^QTG ^^-^{^- VA2-3m) . (11) 

In a singular universe, it is then inevitable that at some instant the (divergent) effective energy 
density violates the entropic threshold of Eqs.(lO) or (11) for a finite value of q. Hence in the 
considered GB gravity (with A > 0) and QTG models of /x < 0, the cosmological singularity 
lies in a region of space-time which is already unphysical, for the apparent horizon entropy is 
negative. 

4 An example of QTG cosmology 

In order to describe the changes in the evolution of the homogeneous and isotropic Universe 
caused by the cubic QTG terms, we next address the problem concerning the construction 
of analytic solutions of the modified Friedmann equations (5) in the particular example of the 
matter stress-energy tensor T^^. , whose components are related by the following linear barotropic 
equation of state [10-13]: 

Po/go = ujoigo) ; tdoigo) = wi-w2/go ■ (12) 

The constant W2, if positive, represents an energy density which contributes to the pressure 
Pq independently of the variable energy density ^0 " thus when the former dominates (i.e. 
^2/ go ^ 1)) the dynamics is driven by a constant energy density with EoS po = —W2, resulting 
in a de Sitter geometry. In this sense, W2 may be identified with the cosmological constant and 



(12) is an example of a simple quintessence model. The constant wi, which may be seen as the 
"matter EoS parameter" (as opposed to the "cosmological constant parameter" ^2) is equal to 
the velocity of sound in the fluid: Vq = dpo/dgo = wi (in Plank units, c= 1 = h). The causality 
condition ti^i < 1 excludes superluminal velocities, and if we also assume that the universe does 
not have a "phantom" phase, then we have to impose \wi\ < 1. 

The dynamics of the universe in such QTG cosmology may be found by solving the differential 
equation for H{t) obtained by combining Eqs.(5) with the EoS (12): 



(1 - 3iiL^H^)H 
the integration of which yields 



Ih\1 - fiL^H^){l + wi) + —W2, 



T 



(13) 



tiH) 



' 3{l+wi) TJJl 



2v 



arctg 



^-«Uarctg(^ 



-Xlog ( gg;g'+g ) + 7^ arccth(i:?/i/A) - 27ri; 



(14) 



Here Hj^ is the only real root of the cubic equation (8) when written in terms of H^ = k^q/Q 
and solved for H^, with Hq = k'^qq/Q considered known. The real and imaginary parts of the 
remaining roots, h^ and its complex conjugate /i^, are denoted by ^ and C,, while x &iid v are 
the real and imaginary parts of a/2 = 2h(h'^~H^)(h^-h^) ' ^o^i*^^ ^^^^ ^^ have chosen one specific 
singular solution H{t — )■ 0) — )■ 00 of Eq.(13), representing big-bang space-times with singularity 
at t = 0. Depending on the initial conditions imposed on H{t) (or equivalently on g{t)) one 
can construct other non-singular bounce- like solutions (both for QTG or EH models), which are 
however out of the scope of the problems discussed in the present paper. 

The function t{H) (14) is not trivially invertible; nevertheless it allows to describe the prop- 
erties of the different periods of acceleration and deceleration of the QTG corrected Universe 
evolution, that can be obtained directly from Eqs.(13), (5b) and (6). We next recall that one 
can also use the deceleration parameter q = —a a/ a? (for a ^ 0), written in a suggestive form: 



q=\{l + ?,p/Q), 



(15) 



in order to determine whether the universe undergoes accelerated {q < 0) or decelerated {q > 0) 
expansion. Thus, the ratio p/ g = cjcfi between the components of the effective energy-momentum 
tensor plays the role of an ^effective equation of state'. Its explicit form : 



l^esiQ) 



-1 + 



(1 + wi) (1 - /iL^x^^V36) g 



W2 



^(1-^^4x4^2/12) 



(16) 



is derived by substituting Eq.(12) into Eqs.(5b) and (6). As expected, for /i = we get Weflf = 
uiQ. Since g ^ iff Weff ^ —1/3, it is convenient to consider only Eq.(16). Observe that as 
one approaches the initial singularity and g diverges, we can take the zeroth order limit of 
l/du-L^x^g"^) ^ 1 in order to show that 1 -|- ujcs ~ ^(1 -|- ujq). Therefore, even if ojq is in the 
"most decelerated range" possible, viz. uq ^ 1, in the considered QTG cosmology we have an 
accelerated phase: Wcff ^ —1/3. Thus, the addition of the cubic QTG terms (1) to the EH 
action results in a new acceleration period at the beginning of the universe. 



'^Changes in the EoS due to the higher curvature terms arising within the context of /(_R)-modified gravity 
have been studied in Refs. [27,28]. 
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Figure 1: Evolution of the equation of state parameter for /i = —10^^: the sohd hne corresponds 
to the effective EoS of QTG as a function of the effective energy density, Eq.(16); the dashed 
black line corresponds to the EH case, a;o(^o)j the red lines are the entropy densities of the 
apparent cosmological horizon, Eq.(3), for (i) A = —9.00 x 10^^ ; (ii) A = —2.50 x 10^^ ; (in) 
A = 7.53 X 10^9. 



This effect may be seen in Fig.l, where ojq{qq)^ given by Eq.(12), is shown as the dashed line, 
while the continuous black line depicts the effective equation of state (16). Here wi = 0, and W2 
is fixed by the observed value of the cosmological constant (cf. Eq.(20)). One can clearly see 
that as the energy density increases the plot of a;cff sinks beneath the line of —1/3, indicating the 
new period of accelerated expansion. For smaller densities, however, there is very little difference 
between the EH plot and the QTG one. Although some of the indispensable features of inflation 
- slow-roll, for example - are not present, we shall freely nominate this early acceleration period 
as an "inflationary" . The duration of such a "rustic inflation" evidently depends on the value 
of /i. By imposing WeffC^/acc) = —1/3, we get a cubic equation: 



(«)l-l)_4 



36 



>c [iL 



t/acc 



^^^acc + «;2 







(17) 



whose positive real roots give the threshold of the two acceleration periods now present in the 
dynamics - the initial one due to QTG and the final due to the cosmological constant. The 
number of distinct real solutions depends on the sign of the discriminant of (17) and, for /i < 
and w\ > —1/3, it determines a critical value 



l/^a 



16(3tt;i + 1)3 
81(1 - wi)x^w^L^ 



for which the discriminant vanishes. Thus if |^| < l/^acd there are two positive real roots for 
(17), corresponding to two distinct periods of acceleration. But for large enough values of 
the gravitational coupling, namely |/i| > |/Xacc|) there is no positive real root to (17) and as a 
consequence the initial QTG acceleration period lasts for such a long time that it merges with 
the final one - the universe is then never decelerated. 

It is worthwhile to remark here that in the QTG counterpart (16) of the EH linear EoS, for 
fx < the effective speed of sound f^g = dp/dg also satisfies the causality conditions — 1 < v'^q < 
1, under the same restrictions on the parameters wi and W2- On the other hand, for /i > 0, we 
see that v'^^ — )• oo as ^ — )• -jj^^f^Jpu. This non-causal behaviour of the effective fluid near the 
points where Co(^) = is yet another reason to consider here only the case of negative values 
of 11. 



Although the form of t(H) given by Eq.(14) does not ahow to analytically determine the 
exact form of the Hubble function H{t), we may invert it in a first order approximation. This 
is possible when the dimensionless quantity \fi\L^H'^ is much smaller than unity. For L = /pi, 
such an approximation is valid during most of the universe history, since H{t) is typically of a 
cosmological order (greater than 1 Mpc ~ 10^^ x /pi). In what follows, we shall refer to this 
approximation as "first order in |/i|": 

H{t;fi) = Ho{t) + fiHi{t) + --- , 
a{t,fi) = ao{t)[l + fiAi{t) + ■■■], 

with the scale factor parametrized as a{t) = e^*^ In zeroth order we have, naturally, the EH 
solution: 

^o(t) = ^^6(rf^Cth [{t - to)/T] , (18) 

2 
3(l+««i) 



oo(t) = aosh'5[(t-to)/T], 5=3^^^. (19) 



Here 1/r = >cy^3{l + wi) W2/8, and oq is a normalization constant. At early times ao{t) ~ 
o-oit/r) the observed behaviour is typical of cosmologies with constant EoS: po/go ~ wi. Later, 
for t/r ^ 1 the cosmological constant 

Ao = }<'^W2/8{l + wi) (20) 

dominates the EoS and the universe enters a final, accelerated, asymptotically de Sitter phase, 
with an asymptotically constant energy density ^a,, = W2/{1 + wi). Then 

ao(t)~exp|2./^(t-to)l • (21) 



In particular, by choosing wi = 0, and thus 6 = 2/3, we see that (12) describes fairly well our 
observed Universe, neglecting inflation and the radiation dominated era: we begin at t = with 
a dust-filled space-time, which ends at a de Sitter space with Aq presenting observed value [10] 

Ao ^ 3.14 X 10^^22 ^ 1-2 ^ (^22) 

if we choose ■W2 accordingly, using Eq.(20). 

The first order corrections can be easily calculated from Eq.(13): 

H,(f) - Hi ft~to _ 1 Y, (^it±ll\ _ 

^^^*) - shH{t-U>)/r)\^^ 8^^[-^) 

-lcth(*^)[sh-^(*^) + i]}; (23) 

A,{t) = 4^{cth (i) [5cth (i) - f ] + sh-^ (i) } , (24) 

where Hi = — x^ 11)2 tL^/ 72(1 + t«i)^. The constant to = 47rf/3(l + wi)fiL^ assures that the 
singularity is placed in t = 0, and as /U — )■ also to — ^ 0. From the scale factor and the Hubble 
function, one can determine all other relevant quantities, and in particular the effective energy 
density Q{t), plotted in Fig. 2. 

As one approaches the initial singularity at t = 0, Eq.(14) shows that H diverges. Eventually 
we then have H > l//pi and the first order approximation is bound to fail - indeed, Hi (t) diverges 
more rapidly than Ho^t) as t — t- 0. This can be clearly seen in Fig. 2: the black continuous line 
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Figure 2: Time evolution of the energy density, for the same parameters of Fig.l. The black 
continuous curve gives the exact function Q{t), obtained from graphical inversion of t(^) given 
by Eq.(14); the black dashed line gives the first order approximation; the dotted gray line 
depicts the EH density go{t); the horizontal gray line marks the final energy density of the early 
acceleration period. 



shows the exact function g{t), obtained from graphical inversion of t{g) given by Eq.(14); the 
black dashed line shows the first order approximation Hq (t) + fiHi (t) . It is clear that the first 
order approximation is only valid valid for times greater than an instant t*, when the curve has 
a maximum, but for t > t* it is in good agreement with the exact solution. To smallest order in 
|/u|, we have 



i*« (2{l + wif/^ + 7r 



V2 



Hi + wi] 



L|^|i/^ 



(25) 



Now, approximating to first order Eq.(17) we find that the initial period ends when the effective 
energy density decreases to the value 



^initial 



l+3wi 7-2 
1-wi 'Pl 



(i^lMl)-^/^ 



which according to Eq.(14) happens at the instant 



U 



(^/3 + 7r)^/2 
3{l + wi) 



L l/i 



1/4 



(26) 



(27) 



This shows that t* is slightly posterior to taco hence an approximation to only first order is not 
sufficient to describe the initial acceleration period created by QTG. 

As we have shown above, if \n\ ^ |/Uacc| we may have an eternally accelerated expansion. Such 
large values of |//|, therefore, do not correspond reasonably to the observed universe. One might 
thus pose the question: What are the restrictions on the values of the gravitational couplings A 
and iJ,, which guarantee the physical consistency of the quasi-topological effects? Regarding the 
initial acceleration period as an infiation era, we may assume that it would occur in the range 
of energies ^pi > £» > 10~^^ x gpi (see [29]), thus we must have /Xacc such that the root of (17) 
lies within this bound. The upper bound of this interval, viz. £»acc ~ 10^^^ x gpi, is indeed 
the case depicted in Figs.l and 2, what serves to demonstrate the validity of the first order 
approximation. Another phenomenological restriction stems from the fact the the apparent 
horizon entropy should not be vanishing for too small values of g. In the GB case, there is no 
initial acceleration period, and this entropic restriction is in fact the only condition we have on 



A. It is quite evident from Eq.(lO) that one may choose A to get ^qb as large as one needs - 
e.g. for A = 3/4 we have ^gb = i?Pi- In QTG instead, as it may be easily seen from Eq.(ll), 
for each given value of fi determining the end of the acceleration period, we can choose A in 
order to place ^qtg: (i) before or (Hi) after ^initial, or even to have (ii) qqtg = initial- This is 
exemplified in Fig.l, where the red lines represent the values of g for which s{q) = 0. Note that 
in cases (ii) and (Hi) the whole initial acceleration period is rendered "unphysical" on account 
of the negative horizon entropy there. 



5 On the late Universe QTG effects 

Although the more significant effects of QTG take place when the curvature, as well as the energy 
density, are big enough - namely, at early times - it turns out that the whole evolution of the 
universe is modified by the QTG terms. Clearly, at later times, as the curvature diminishes, 
the cubic and quadratic terms in (1) become more and more negligible and the first order 
approximation made in the last section is then justified. 

An example of such changes is given by the fact that the cosmological constant, which 
characterizes the geometry of the asymptotically dS4 spaces in the limit t — )■ oo, is not equal 
to the bare one, Aq, defined by the matter Lagrangian. Indeed, in QTG one observes the 
effective cosmological constant A^s, related to the Hubble function H{t). To first order in ^, 
one may determine Aeg by the approximation (23) for H{t), or else by inverting directly the 
exact equation 



Ao = Aeff (1 - fiL^Als/9) , 



(28) 



obtained from Eq.(8): Acs ~ Aq + -L^//Aq/9 . This is in fact a general result, valid for every 
asymptotically de Sitter space, regardless of the particular matter EoS leading to the final dS4 
vacuum. Notice that for the very small value of the observed cosmological constant, both Aq 
and Acfi are practically equal. 

Another feature of asymptotically de Sitter space-times is the presence of a future event 
horizon. Its comoving radius is given by the integral rj(t) = J^ dt/a{t) and may be calculated 
to first order in fi by using the results (19) and (24). At zeroth order (i.e. in the EH case) this 
yields 



rf\t) 



Sao ch* (t/r) ^ ^ 

while the first order QTG correction is given by: 
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5 1+5. 2+5. 1 
2' 2 ' 2 ' ch^ (i/r) 



(30) 



The function rf{t) describes the past light cone of the (infinite) future of the comoving observer 
at the origin. The tip of this cone is placed at {rj = 0, t = 00}. Then as t — t- 00 the physical 
radius, lf{t) = a{t)rf{t) becomes equal to the (constant) de Sitter radius 1/H. The first order 



approximation to lf{t) can be easily obtained with the help of Eq.(24): for example its value 
today, at t = 13.8 Gyrs, is If ^ (4.90 x 10^ + ^ x 9.10 x 10~248) Mpc . 

Alternatively, the comoving radius rp(t) = J^ dt/a{t) describes the future light cone for an 
observer at the origin, starting from its tip at {r^ = 0, t = t}. In a singular universe, if this 
tip is placed at the beginning of time, t = 0, then rp{t) is the particle horizon. In practice, the 
choice of t determines only a constant, for we may write 



rp{t) = - / dt/a{t) + / dt/a = -rf{t) + r/(t) . (31) 

Jt Ji 

Thus for example, in the EH case with t = 0, Eq.(29) gives 

rp li; - rhf^ftM^^^ 2' 2 ' 2 ' ch^(t/T) +2ao^l2' 2 J' ^'^^) 



ch" (t/r) 



ch^ (t/r) 



and in general Eqs.(29) and (30) determine also the first order correction for rp{t), for a given t. 

Recall that the first order approximation for the scale factor is only valid for t > t^,, with 
tif denoting the instant where the approximation fails, Eq.(25). Therefore in QTG we have the 
"technical" impossibility of placing the tip of the light cone on the initial singularity - the earlier 
we may place it is at t = t*. Due to the fact that t^, (at our first approximation) is slightly 
posterior to taco we are then technically prevented from describing the increase of the particle 
horizon during this early "inflationary" period. Consequently, we cannot determine its number 
of e-foldings, and whether it solves the usual problems of non-inflationary cosmology - such as, 
for example, the horizon problem. 

The knowledge of the above expressions for the radii rp(t) and rf{t) provides an analytic 
description of the causal diamond [30,31] and its comoving volume K;.d.(Oi ^^ the QTG cosmo- 
logical model under investigation. The former is defined as the intersection of the causal past 
of the "point" {rj = 0, t = oo} and the causal future of {vp = 0, t = t}. At each instant t the 
volume is given by Vcd. = ^TTr^{t) with r = r^ on the upper and r = rj on the lower half of 
it. It is easily seen that I^.d. has a single maximum at t = tedgej when rp(tcdgc) = ^/(^edge) [31]. 
With the aid of Eq.(31), the instant tcdgc is determined from the equation rj(tcdgc) = 2^f{t)^ 
which is exact, i.e. independent from the first order approximation. If rf[t) is a continuous 
function, then we conclude that a first order correction to t must imply a first order change in 
tedge- Therefore for small t* = t in QTG, the edge of the causal diamond occurs at an instant 
diff'ering not more than at first order from the corresponding value in EH. 

In EH gravity with linear EoS (12) and wi = 0, describing an early universe dominated by 
dust®, Bousso et al. [10] have calculated V^.d.(^) ^"^^ used it to predict the order of magnitude of 
the observed cosmological constant. Their analysis, based on the rather universal (phenomeno- 
logical) evaluation of the entropy production rate in our universe, demonstrates that only when 
Aq ~ 10^^^^ X /pj the causal entropic principle (CEP), requiring maximal entropic production 
within the corresponding causal diamond volume, is fulfilled. According to their arguments the 
main production of bulk entropy occurs during the matter-dominated era. Therefore one can 
perform a similar analysis for the QTG extension of this linear EoS model, by imposing the 
CEP conditions for the causal diamond whose inferior tip is placed at t = t*, thus respecting 
the restrictions of the considered first approximation. Surely, even if the tips of the causal di- 
amonds in QTG and in EH gravity were placed at the same point - say, by replacing t = 
with t = t^, in the EH case as well -, still they would not be identical due to the changes in 



*In which not only the initial inflationary period, but also the radiation-dominated era of the concordance 
model are absent. 
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the dynamics given in (30). However, as we saw above, the first order QTG corrections we 
are considering do not change significantly the comoving volume of the causal diamond, and in 
fact they lead to the same prediction for the magnitude of the effective cosmological constant 
Aeff « Ao + LV Ao/9 ~ 10-122 X /pf. 

Throughout this discussion, we have not taken into account the entropic bounds derived in 
Sect. 3. Regardless of our technical restrictions for placing the inferior tip of the causal diamond, 
in the case of ^ < 0, we cannot place it before the instant tqTG when the apparent horizon's 
entropy vanishes. The same is true in the GB case, for A > 0. We are however assuming that 
iqTG is very small, in particular that tqxG ^ i*- This can always be set for an appropriate value 
of A. Either way, the entropic restrictions do not allow the inferior tip of the causal diamond to 
be placed on the singularity. 

Let us list in conclusion a few open problems concerning the considered (linear EoS) cos- 
mological model of the cubic Quasi- Topological gravity: (a) The stability conditions for the 
asymptotic dS4 cosmological QTG solutions, which requires the calculation of the spectrum of 
the corresponding linear fiuctuations at least in the probe approximation; (b) The analysis of 
the properties of QTG models with more realistic matter content by considering EoS or equiv- 
alently (string inspired) matter superpotentials [32] giving rise to early-time infiation with a 
desired slow-roll behaviour, which is under investigation [14]. It is worthwhile to also mention 
that the methods and some of the results of the present paper seem to have a straightforward 
"cosmological" application to the recently constructed "higher curvature" quartic QTG exten- 
sion [33] of the EH gravity, as well as to the case of spatially curved FRW solutions of the 
considered four dimensional, cubic QTG model. 

Acknowledgments. We are grateful to C.P. Constantinidis for his collaboration in the initial 
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